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Abstract. A variant of the quantum marginal problem was known from early sixties 
as N-representability problem. In 1995 it was designated by National Research Council 
of USA as one of ten most prominent research challenges in quantum chemistry. In 
spite of this recognition the progress was very slow, until a couple of years ago the 
problem came into focus again, now in framework of quantum information theory. In 
the paper I give an account of the recent development. 
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1. Introduction 

The quantum marginal problem is about relations between spectrum of mixed state pab 
of two (or multi) component system TCab = 'Ha ®'Hb and that of reduced states pA 
and pb of the components. This problem emerges a couple of years ago in framework 
of quantum information theory. 

The problem can be stated in plain language as follows. Let M = [rriijk] be complex 
cubic matrix and Mi, M2, M3 be Gram matrices formed by Hermitian dot products of 
parallel slices of M. We are seeking for relations between spectra of these matrices. 

In section 3 we survey some recent results that laid the ground of the quantum 
marginal problem and discuss a solution of this problem based on geometric invariant 
theory. 

A variant of QMP dealing with system of fermions, like electrons in an atom or a 
molecule, was known from early sixties as A^-representability problem (Coulson, 1960). 
Its solution allows to calculate nearly all properties of matter which are of interest for 
chemists and physicists (Coleman and Youkalow 2000). By this reason it was designated 
as one of ten most prominent research challenges in theoretical chemistry (Stillinger et 
al 1995). In section 3 we outline a general solution of the A^-represent ability problem 
for one particle reduced density matrix. For systems of rank < 8 explicit sets of linear 
constraints are given for pure state A^ representability. A representation theoretical 
interpretation of A^-representability plays crucial role in the calculations. 

The reduction pab ^ Pa is known in mathematics as contraction. For example, 
Ricci curvature of a Riemann manifold is a contraction of its Rienmann curvature. 
The results of this paper imply some inequalities between spectra of Riemann and 
Ricci curvatures, see Example 14.2.31 Recall that in general relativity Ricci curvature is 
governed by the energy-momentum tensor, i.e. by physical content of the space, while 
Riemann curvature is responsible for its geometry and topology. The above constraints 
impose some bounds on influence of matter on geometry. 

2. Classical marginal problem 
2.1. Marginal disributions 

Let's start with classical marginal problem which asks for existence of a "body" in M"' 
with given projections onto some coordinate subspaces C M", / C {1, 2, . . . , n}, i.e. 
existence of probability density p{x) = p{xi,X2, ■ ■ ■ , Xn) with given marginal distributions 

Pi{xi) = / p{x)dxj,J = {l,2,...,n}\I. 

In discrete version the classical MP amounts to calculation of an image of a 
multidimensional symplex, say A = {pijk > 0| YliPijk = 1}) under a linear map like 

Pijk ^ {Pij-,Pjk-,Pki)i 
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Pij / ^ Pijkj Pjk / j Pijki Pki / _, Pijk- 

k i j 

The image 7r(A) is a convex hull of the projections of vertices of A. So the classical 
MP amounts to calculation of facets of a convex hull. In high dimensions it might be a 
computational nightmare (Pitowsky 1989). 



2.2. Classical realism 

Let Xi : Ha "Ha be observables of quantum system A. Actual measurement of 
Xi produces random quantity Xi with values in Spec (Xj) and density Pi{xi) implicitly 
determined by expectations 

u{x,)) = m{xm 

for all functions / on spectrum Spec(Xj). For commuting observables G / the 
random variables Xi,i G / have joint distribution pi{xi) defined by similar equation 

{f{xi)) = m{Xjm, V/. (2.1) 

Classical realism postulates existence of a hidden joint distribution of all variables Xi. 
This amounts to compatibility of the marginal distributions 1)2.11) for commuting sets 
of observables X/. Bell inequalities, designed to test classical realism, stem from the 
classical marginal problem. 

2.2.1 Example. Observations of disjoint components of composite system TY^ ® T^b 
always commute. For two qubits with two measurements per site their compatibility is 
given by 16 inequalities obtained from Clauser-Horne-Shimony-Holtmequality (Calauser 
et al 1969) 

(ai6i) + (0261) + (0262) - (ai&2) + 2 > 

by spin flips \—>- ±aj and permutation of the components A ^ B. Here {aibj) is 
expectation of product of spin projections onto directions i, j in sites A,B. 

2.2.2 Example. For three qubits with two measurements per site the marginal constraints 
amounts to 53856 independent inequalities (Pitowsky et al 2001). This example may 
help to disabuse us from overoptimistic expectations for the quantum marginal problem 
to be discussed below. 

2.2.3 Example. Univariant marginal distributions are always compatible, e.g. we can 
consider Xi as independent variables. However under additional constraints, say for 
a "body" of constant density, even univariant marginal problem becomes nontrivial. 
For its disctete version Cale-Ryser theorem (Gale 1957) tells that partitions A,/i are 
margins of a rectangular 0/1 matrix iff majorization inequality A -< /x* holds. Here 
marginal values arranged in decreasing order are treated as Young diagrams 



A = (5,4,2,1) = 



A* = (4, 3, 2, 2,1) = 
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/X* stands for transposed diagram, and the majorization order X ~< u is defined by 
inequalities 

Ai < z/i 

Ai + A2 < 1^1 + V2 

Ai + A2 + A3 < vi + 1^2 + 1^3 



3. Quantum marginal problem 

3.1. Reduced states 

Density matrix of composite system AB can be written as a linear combination of 
separable states 

PAB = ^ao,PA(^ P%, (3.1) 

a 

where p% p% are mixed states of the components A, B respectively, and the coefficients 
Cq, are not necessarily positive. Its reduced matrices or marginal states may be defined 
by equations 

PA = Y,a^Ti{p^B)pl :=TrB(pAB), 

a 

PB = $^a,Tr(p>^ ■.= 1:ia{pab). 

a 

The reduced states Pa, Ps are independent of the decomposition (|3.ip and can be 
characterized intrinsically by the following property 

{XA),,,=TripABXA)=TTipAXA) = {XA),,, V Xa, (3.2) 

which tells that p^ is a "visible" state of subsystem A. This justifies the terminology. 
3.1.1 Example. Let's identify pure state of two component system 

with its matrix [ipij] in orthonormal bases ai,l3j of Ha,'Hb- Then the reduced states of 
ip in respective bases are given by matrices 

Pa = i^H, pB = 'ip'ip^ (3.3) 

which have the same non negative spectra 

Specp^i = Specps = A (3.4) 

except extra zeros if dimTCA 7^ dimTi^- The isospectrality implies so called Schmidt 
decomposition 



^ = $^V^^.^®^f, (3.5) 

i 

where ipf, ipf are eigenvectors of Pa, Pb with the same eigenvalue Aj. 
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In striking contrast to classical case marginals of a pure state ^ iI'a^4'b are mixed 
ones, i.e. as Srodinger 1935 puts it "maximal knowledge of the whole does not necessarily 
includes the maximal knowledge of its parts. " He coined the term entanglement just to 
describe this phenomenon. 

3.2. Statement of the problem 

Quantum analogue of the classical marginal distribution is reduced state pA of composite 
system Hab = 'Ha ^Ti-B- Accordingly, most general quantum marginal problem asks 
about existence of mixed state pi of composite system 

with given reduced states pj for some J C I (cf. with classical settings section 
Additional constraints on state pj may be relevant. Here we consider only two variations: 

• Pure quantum marginal problem 

dealing with marginals of pure state pj = and more general 

• Mixed quantum marginal problem 

corresponding to a state with given spectrum A/ = Specp/. 

Both versions are nontrivial even for univariant margins (cf. Example 12. 2. 3^ . In 
this case reduced states pi can be diagonalized by local unitary transformations and 
their compatibility depends only on spectra Aj = Specpj. Note that mixed QMP say 
for two component system Has = T^a ® ?^_b is formally equivalent to the pure one for 
three component system Hab ® Ha ®Hb- 

Pure quantum marginal problem has no classical analogue, since projection of a 
point is a point. For two component system Ha ®Hb marginal constraints amounts to 
isospectrality Spec p^i = Specps, see Equation ()3.4j) . For three component system the 
problem can be stated in plain language as follows. 

Problem 3.2.1. Let A = [Aijk] be complex cubic matrix and Ai,A2,As be Gram 
matrices formed by Hermitian dot products of parallel slices of A. The question is 
what are relations between spectra of matrices Ai, A2, A3? 

Unfortunately methods of this paper can't be applied directly to overlapping 
marginals like pab, Pbc, Pca- 

3.3. Some known results 

Here are some recent results that laid the ground of the quantum marginal problem. 
They all emerge from quantum information settings in two or three years. 

Theorem (Higuchi et al 2003). For array of qubits <S)'i=i'H.i, dim?-^j = 2 all 

constraints on margins pi of a pure state are given by polygonal inequalities 

A. < E 
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for minimal eigenvalues Aj of pi. 

This result was proved independently by Sergey Bravyi who also managed to crack 
mixed two qubit problem. 

Theorem (Bravyi 2004). For two cubits Ha'Si'Hb solution of the mixed QMP is given 
by inequalities 

min(A^,AB) > A^^ + A^^, 

A^ + Xb > + A^^ + 2Xf^ 

\Xa - Xb\ < mm{Xt^ - A^^, A^^ - Af^), 
where Xa,Xb are minimal eigenvalues of pA, Pb and X^^ > X^^ > Xf^ > Xf^ is 
spectrum of pab- 

Finally for three qutrits the problem was solved by Matthias Pranz using rather 
advanced mathematical technology and help of a computer. An elementary solution was 
found independently by Astashi Higuchi. 

Theorem (Franz 2002, Higuchi 2003). All constraints on margins of a pure state 
of three qutrit system Ha ^Hb® Hq are given by the following inequalities 

A^ + A? <A^ + A? + A^ + AJ, 

Xl + Xl < A^ + A^ + A^ + AJ, 

^^3 + ^^2 < A2 + AJ + A3 + A2, 

2X1 + A? < 2X\ + Aj + 2X1 + 
2X1 + A2 < 2X\ + A? + 2X1 + ^2, 

2A« + A" < 2X\ + A^ + 2X1 + ^2, 
where a, 6, c is a permutation of A,B, C , and the marginal spectra are arranged in 
increasing order Ai < A2 < A3. 

Note that in contrast to the classical marginal problem, linearity of the quantum 
marginal constraints is a surprising nontrivial fact. 

3.4- Main theorem 

A general solution of the quantum marginal problem, based on geometric invariant 
theory, has been found recently (Klyachko 2004). 

Theorem 3.4.1. For two component system H-ab — Ha <8) Hb of format m x n all 
constraints on spectra X^^ — SpecpAB, X^ — Specp^; A^ = Specpe arranged in 
decreasing order are given by linear inequalities 

i j k 

where a : ai > 0,2 > ■ ■ ■ > am, 6 : 61 > 62 > ■ ■ ■ > "^CLi = = are "test 

spectra", a + 6 = {oj + bj}l, and u G S„i, v G Sn, w G S„m are permutations, subject to 
crucial condition c^^(a, 6) 7^ 0/ topological nature to be explained below. 
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3.4-2 Remark. Here {a + b)k denotes k-th term of the sequence + bj arranged in 
decreasing order. The coefficient c^„(a, b) depends only on the order in which quantities 
+ bj appear in the spectrum (a + The order changes when the pair (a, b) crosses 
hyperplane 

Hij\ki '■ <Xi + bj = ak + bi. 

The hyperplanes cut the set of all pairs (a, b) into finite number of pieces called cubicles. 
For each cubicle one have to check inequality ()3.6|) only for its extremal edges. Hence the 
marginal constraints amounts to a finite system of inequalities, but the total number of 
extremal edges increases rapidly. Here are some sample data for arrays of qubits. 



# qubits 


2 


3 


4 


5 


6 


# edges 


2 


4 


12 


125 


11344 



Unfortunately for most systems the marginal constraints are too numerous to be 
reproduced here. Therefor we give only summary table, which shows how complicate 
may be the answer. 



System 


Rank 


7^ Inequalities 


2x2 


2 


7 


2x2x2 


3 


40 


2x3 


3 


41 


2x4 


4 


234 


3x3 


4 


387 


2x2x3 


4 


442 


2x2x2x2 


4 


805 



3.5. Hidden geometry and topology 

Here we explain the meaning of the coefficient c^^{a, b) in statement of the theorem and 
show how they can be calculated. Consider set of Hermitian operators Xa : Ha 'Ha 
with given spectrum Spec(X^) = a and call it flag variety 

^aiH-A) ■= {Xa I Spec(X^) = a}. 
For two flag varieties J^aiH-A) and J^biHs) define map 

fab '■ ^a{HA) X J^biHs) ^ Ta+bO^A ® Hb), 

XaX Xb I — ^ Xa01 + 1®Xb. 
The coefficients c^„(a, 6) come from the induced morphism of cohomology 

fl, : H\:Fa+b{HAB)) ^ H*{j^a{nA)) ® H*{j^b{nB)) 

written in the basis of Schubert cocycles 
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We'll give below an algorithm for their calculation. For this we need another description 
of the cohomology of flag varieties due to (Bernstein et al 1973). Specifically, for simple 
spectrum a eigenspaces of operator Xa G J-'aiTi-A) of given eigenvalue Oj form a line 
bundle Cf on J^aO^A)- Their Chern classes xf = Ci{Cf) generate cohomology ring 
H*{Ta{'HA)) and in this setting morphism (yj*^, admits a simple description 

where (a + h)k = ctj + hj. In terms of the canonical generators Xi = Ci(£j) Schubert 
cocycle is given by Schubert polynomial (Macdonald 1991) 



n—ln—2 

X n 



■ X. 



n-1), 



where w G S'„ is a permutation of 1, 2, . . . , n, wq = {n,n — 1, ... ,2,1) reversion of the 
order, operator d^, = di^di^ ' ' '^in is defined via minimal decomposition w = Si-^Si^ ' ' ' 
into product of transpositions Si = {i,i + 1), ^ = i{w) is the number of inversion of w 
called its length , and finally 



/(• • • Xi^i, . . .) /(..., Xi-^-l, Xi, . . .) 



is finite difference operator. This leads to the computational formula 



cZia,b) = d^d^SUx 



„AB\ 



(3.7) 



where (a + b)k = + bj and (i{w) = i{u) + i{v). It can be easily implemented into a 
computer progrm. Recall that in order to get a finite system of inequalities one have 
also to find all the extremal edges and use them as the test spectra (a, b). 

3.5.1 Example. Note that for identical permutations u,v,w the coefficient c^^(a,6) is 
equal to 1. Hence inequality 

Y^aX + Y,b,\f <Y,{a + b),Xr 

i j k 

holds for all test spectra (a, b). This amounts to finite system of basic inequalities (Han 
et al 2004) 

+ + + <\^'' + \r + --- + ^kn. k<m. 



Af + Af + 



Af < A^^ + + h A;^f , 



< n. 



3. 6. Array of qubits 

The calculations can be essentially reduced using the following result. 

Theorem 3.6.1. In the setting of Theorem \3.4-l\ all marginal constraints are given by 
inequalities \3. 6|) with c'^^{a,b) = 1. 

3.6.2 Example. One can check that for array of qubits ^^^I'Hi all inequalities ()3.6p 
with odd coefficient c^^„2...„„(ai, 02, • • • , o-n) can be obtained from the basic inequality 

^ ai(Af^ - A^'^) < ^(±ai ± as ± ■ ■ ■ ± a„)fcAfc 
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by transposition A^+i, k = odd, combined with sign change of a summand in 

LHS. Here aj > — are the test spectra, A and A*^*-* are spectra of state p of the system 
and its one qubit reduced state p*^*^ respectively. To get a finite system one has only to 
find the extremal edges. For large this may be a challenge, see Remark 13.4.21 

3. 6. 3 Example. For 3-qubit the theorem returns the following list of marginal inequalities 
grouped by their extremal edges. The first inequality in each group is the basic one. 
The transposed eigenvalues in modified inequahties typeset in bold face. We expect 
Aj = A['^ — A2 ^ to be arranged in increasing order Ai < A2 < A3. 







A3 


< 


Ai + A2 + A3 + A4 — A5 — 


■Ae- 


-A7- 


As. 




A2 


+ A3 


< 


'2\\ -\- 2A2 — 2A7 — 2A8. 










+ A2 


+ A3 


< 


3Ai + A2 + A3 + A4 — A5 


-Ae 


-A7- 


- 3A8, 


-Ai 


+ A2 


+ A3 


< 


3A2 + Ai + A3 + A4 - As 


-Ag 


i — A7 


-3A8 


-Ai 


+ A2 


+ A3 


< 


3Ai + A2 + A3 + A4 - As 


-Ae 


-As 


-3A7 


Ai- 


h AaH 


h2A3 


< 


4Ai + 2A2 + 2A3 - 2A6 - 


2A7 • 


-4A8, 




-Ai- 


h AaH 


h2A3 


< 


4A2 + 2Ai + 2A3 - 2A6 - 


-2A7 


-4A8 




-Ai - 


h AaH 


h2A3 


< 


4Ai + 2A2 + 2A4 - 2A6 - 


2A7 


-4A8, 




-Ai- 


h AaH 


h2A3 


< 


4Ai + 2A2 + 2A3 - 2A5 - 


2A7 


-4A8, 




-Ai- 


h AaH 


h2A3 


< 


4Ai + 2A2 + 2A3 - 2A6 - 


2A8 


-4A7 





4. N-representability problem 

Ji^.l. Physical background 

The quantum marginal problem may be complicated by additional constraints on state 
ip. For example, Pauli principle implies that state space of identical particles shrinks 
to symmetric tensors 5*^7^ C T-i®^ for bosons and to skew symmetric tensors A^Ti for 
fermions. For such systems reduced density matrices appear in the second quantization 
formalism in the form 

p(i) = {ip\a\aj\ilj) = 1 particle RDM, 
p(^) = {iplalojakailip) = 2 particle RDM, etc. 

Their physical importance stems from the observation that, say for fermionic system, 
like multi electron atom or molecule, with pairwise interaction 

N 

H = ^Hi + Y^ Hij 

i i<j 

the energy of state ip depends only on 2-point RDM 
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where if*^^^ = j^[Hi + H2\ + H12 is reduced two particle Hamiltonian. This allows, for 
example, to express the energy of ground state £^0 via 2-point RDM 




The problem however is that it is not obvious what conditions the RDM itself should 
satisfy. This is what the quantum marginal problem is about. In this settings it is known 
from early sixties as N- representability problem (Coulson 1960, Coleman 1963). Later 
in mid 90-th the problem was regarded as one of ten most prominent research challenges 
in quantum chemistry (Stillinger et al 1995). Its solution allows to calculate nearly all 
properties of matter which are of interest to chemists and physicists. For current state 
of affairs and more history see (Coleman and Yukalov 2000, Coleman 2001). 

4.2. One point N -representability 

Here we outhne a solution of the problem for one point reduced states. Following 
chemists we treat them as electron density and accordingly use normalization Tr p^^^ — 
N while keeping Trp = 1. There are few cases where complete solution of one point 
A'"-representability was known prior 2005: 

• Pauli principle: < Aj < 1, A = Spec p*^^^. This condition provides a criterion 
for mixed N -representability (Coleman 1963). 

• Criterion for pure A'"-representability for two particles A^Hr or two holes A^^'^Hr 
is given by even degeneration of all eigenvalues of p'^^\ except (resp. 1) for odd 
r = dimHr (Coleman 1963). 

• For system of three fermions of dimension six A^Tie all constraints on one point 
reduced matrix of a pure state are given by the following (in) equalities 



where Ai > A2 > A3 > A4 > A5 > Ae is spectrum of p^^\ 

The last result belongs to Borland and Dennis 1972 who commented it as follows: 

"We have no apology for consideration of such a special case. The general 
N -representability problem is so difficult and yet so fundamental for many 
branches of science that each concrete result is useful in shedding light on the 
nature of general solution. " 

For more then 30 years passed after this theorem no other solution of A^-representability 
problem has been found. Borland and Dennis derived their criterion from an extensive 
computer experiment, and later proved it with help provided by M.B. Ruskai and 
R.L. Kingsly. They also conjectured solutions for systems A^Ht, A^TYy, A^Tig, e.g. for 



Al + Ag — A2 + A5 — A3 + A4 — 1, A4 < A5 + Xq, 



(4.1) 



A^Hr one point pure representability is given by 4 inequalities 

Ai+A6 + A7>l, A2 + A5 + A7>1, 
A3 + A4 + A7>1, A3 + A5 + A6>1, 



(4.2) 
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but they failed to prove them. The conjectures turn out to be true and covered by the 
following general result. 

Theorem 4.2.1. For mixed state p of n -fermion system A"7ir (^H constraints on spectra 
V = Spec p and X = Spec p^^^ are given by inequalities 

^aiA^(i) < (4-3) 
« j 

for all "test spectra" a : Oi > 02 > ■ ■ ■ > ctr? ^ = 0. Here A'^a = {aj^+aj2 + - ■ ■+aj„} | 
consists of all sums + + ■ ■ ■ + a^^, zi < Z2 < ■ ■ ■ < arranged in decreasing order, 
and V G Sr,w G S^n^ are permutations subject to topological condition c^{a) ^ to be 
explained below. 

4-2.2 Remark. Recall that the spectra A and u are arranged in decreasing order and 
normalized to trace n and 1 respectively. Similarly to Theorem 13.4.11 the coefficients 
c^(a) are defined via flag variety J-'a{'H) := {A : H ^ H \ Spec(A) = a} and morphism 

where operator A^"^ : A"7^ — > A"7^ acts as differential 



A^"^ : Xi A 0:2 A . . . A a;„ ^ Xi A X2 A . . . A Axi A 



A Xr 



The coefficients c^(a) come from the induced morphism of cohomology 

ifl : ff*(J-^.,(A"7^)) ^ H*{Ta{H)) 
written in the basis of Schubert cocycles 

They can be calculated by equation 

where + + ■ ■ ■ + aj„ is k-th term of the sequence A'^a, cf. section 13.51 

4.2.3 Example. For system A^7i4 the marginal constraints on u = Spec p and A = 

Spec p^^^ are given by inequalities 

2Ai < 1^1 + 1^2 + 1^3 
2A4 > 1^4 + 1^5 + 1^6 
2(Ai - A4) < z/i + z/2 - - t'e 
Ai + A2 - A3 - A4 < z^i - z/6 (4.4) 
Ai - A2 + A3 - A4 < min(z/i - z/5, z/2 - z/e) 
|Ai - A2 - A3 + A4I < min(z/i - z/4, z/2 - z/5, z/3 - z/e) 
2 max(Ai - A3, A2 - A4) < min(z/i + z/3 - z/5 - z/g, z/i + z/2 - z/4 - z/g) 
2 max(Ai - A2, A3 - A4) < min(z/i + z/3 - z/4 - z/g, z/2 + z/3 - z/5 - z/g, z/i + z/2 - z/4 - z/5). 
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In these inequalities we use standard normalization Tr p = Tr p^^^ . Recall that all spectra 
are written in decreasing order. 

4-2.4 Example. Similar compatibility conditions for system A^Tis contain 460 
independent inequalities which can't be reproduced here. 

4.2.5 Remark. In tensor algebra the reduction pab ^ Pa is known as contraction. 
Most mathematicians are familiar with this procedure from differential geometry, where 
say Ricci curvature tensor Ric : T — > T is defined as contraction of Riemann curvature 
R : A^T — A^T (here T stands for tangent bundle). This allows to interpret inequalities 
fl4.4|) as relations between spectra of Riemann and Ricci curvatures of four-manifold. 
Recall that in general relativity Ricci curvature is governed by energy-momentum tensor, 
i.e. by physical content of the space, while Riemann curvature is responsible for its 
geometry and topology. The above constraints impose some bounds on influence of 
matter on geometry. 

4.3. Pure N -representability in dimension < 8 

Here I'll give an account of our joint work with Murat Altunbudak. The details will be 
pubhshed elsewhere. 

Formally solution of pure marginal problem can be deduced from inequalities ()4.3p 
of Theorem 14 . 2 . 1 1 by putting z/j = for z 7^ 1. However for a system like A^Tig we confront 
with an immense symmetric group of degree (^) = 70. A representation theoretical 
interpretation of iV-representability discussed below allows to circumvent this difficulty. 

Let's start with decomposition of a symmetric power of A"7i, called plethysm, into 
irreducible components 

S^'^(A-n) =J2^>^'^>^ (4.5) 

A 

of the unitary group UiH). The components Hx, entering into the decomposition with 
some multiplicities > 0, can be parameterized by Young diagrams 

A : Ai > A2 > ■ • • > A, > 

of size |A| = Yli \ = n ■ m that fit into r x m rectangular, r = dimTi. It is instructive 
to treat the diagrams as spectra. We are interested in asymptotic of these spectra as 
m — >■ 00 and therefor normalize them to a fixed size A = A/m, Tr A = n. 

Theorem 4.3.1. Every A obtained from irreducible component Tix C S'^{A"'T-[) is 
spectrum one point reduced matrix p^^^ of a pure state ip G A^Ti. Moreover every one 
point reduced spectrum is a convex combination of the spectra A with bounded m < M. 

A similar result holds in standard settings of the quantum marginal problem (Franz 
2002, Christandl and Mitchison 2004, Klyachko 2004, Christandl et al 2005). 
Note that representation S"^{A'''~^H) is dual to S'™(A"7f) and hence 

5™(A'-"7^) = ^mA7^A*, (4.6) 

A 
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where A* is complement of the diagram A to the rectangle r x m and the multiplicity 
mx is the same as in ()4.5|) . Thus we arrived at the following particle-hole duality. 

Corollary 4.3.2. Marginal constraints on spectrum of reduced matrix of a pure state 
for system A'^^"'7ir can be obtained from that of the system A^Tij. by substitution 
Aj I— >■ 1 — Ar+i-i- 

4-3.3 Example. There are few cases where decomposition ()4.5|) is explicitly known, for 
example 

A|=2m, A=even 

where the sum is extended over diagrams A C r x m with even multiplicity of every 
nonzero row (Macdonald 1995). Together with theorem 14.3.11 and the particle- hole 
duality this implies Coleman's criteria of pure N-representability for systems of two 
particles A^TCr and two holes A^~'^TCr mentioned at the beginning of section 14.21 

4.3.4 Example. Borland-Dennis equations ()4.1|) imply that every component Tix C 
5""(A^7i6) is selfdual A = A*. It seems mathematicians missed this fact, which holds 
only for this specific system. Note that wedge product ensure selfduality of A^Hq and 
hence of the plethysm S'^{A^T-Iq). However apparently there is no simple way to extend 
this to every component Hx C S"^{A^T-Cq). 

Theorem 14 . 3 . 1 1 for any fixed M gives an inner approximation to the set of all possible 
reduced spectra, while any set of inequalities ()4.3p of theorem 14.2.11 amounts to its 
outer approximation. This suggests the following approach to pure A^-representability 
problem, which combines both theorems. 

• Find all irreducible components Tix C 5"'"(A"7i) for m < M starting with M = 1. 

• Calculate convex hull of the corresponding reduced spectra A. 

• Check whether or not all inequality defining facets of the convex hull fit into the 
form ()4.3j) of Theorem 14.2.11 

• If they do then all inequalities are found. Otherwise increase M M + 1. 

4.3.5 Remark. The success of this approach depends on the degrees of generators of the 
module of covariants of the system A^'TCr- Generically the degrees are expected to be 
huge as well as the whole number of the resulting inequalities. However for systems of 
rank r < 8 and for r = 9, n ^ 4,5 the module of covariants is free (Vinberg and Popov 
1992) and the degrees of the generators should be reasonably small. 

Indeed an inexpensive PC, assisted with some dirty tricks, managed to resolve A^- 
represent ability problem for rank r < 8. Recall that for two fermions or two holes the 
answer is known, see section Together with particle-hole duality this bound us to 
the range 3 < n < r/2. The corresponding constraints are listed below. They are 
grouped by the extremal edges and use chemical normalization Aj = n. 

• A^He. 

Ai + Aq = A2 + A5 = A3 + A4 = 1, A4 < A5 + Ag 
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-4Ai + 3A2 + 3A3 + 3A4 + 3A5 - 4A6 - 4A7 < 2 
3Ai - 4A2 + 3A3 + 3A4 - 4A5 + 3A6 - 4A7 < 2 
3Ai + 3A2 - 4A3 - 4A4 - 3A5 + 3A6 - 4A7 < 2 
3Ai + 3A2 - 4A3 + 3A4 - 4A5 - 4A6 + 3A7 < 2 

3Ai — A2 — A3 — A4 — Ag — As — A7 + 3A8 < 1 

— Ai + A2 + A3 + A4 + A5 — Ag — A7 — Ag < 1 
Ai + A2 — A3 — A4 + A5 + Ag — A7 — Ag < 1 
Ai + A2 — A3 + A4 — A5 — Ag + A7 — Ag < 1 
Ai — A2 + A3 + A4 — A5 + Ag — A7 — Ag < 1 

2Ai + A2 - 2A3 - A4 - Ag + Ag < 1 
2Ai — A2 — A4 + Ag — 2A7 + Ag ^ 1 
A3 + 2A4 - 2A5 - Ag - A7 + Ag < 1 
Ai + 2A2 - 2A3 - A5 - Ag + Ag < 1 
2Ai — A2 + A4 — 2A5 — Ag + Ag < 1 

5Ai + 5A2 - 7A3 - 3A4 - 3A5 + Ag + A7 + Ag < 3 
5Ai - 3A2 - 3A3 + A4 + As + SAg - 7A7 + Ag < 3 

5Ai + A2 - 3A3 + A4 - 3A5 + Ag - 3A7 + Ag < 3 
Ai + A2 + A3 + 5A4 - 3A5 - SAg - 3A7 + Ag < 3 
Ai + 5A2 - 3A3 + A4 + As - 3Ag - 3A7 + Ag < 3 

9Ai + A2 - 7A3 - 7A4 - 7A5 + Ag + A7 + 9Ag < 3 
9Ai - 7A2 - 7A3 + A4 + As + Ag - 7A7 + 9Ag < 3 

7Ai - A2 - A3 - A4 - A5 + TAg - 9A7 - Ag < 5 
7Ai - A2 - A3 + 7A4 - 9As - Ag - A7 - Ag < 5 
7Ai + 7A2 - 9A3 - A4 - As - Ag - A7 - Ag < 5 
-Ai - A2 + 7A3 + 7A4 - As - Ag - 9A7 - Ag < 5 
-Ai + 7A2 - A3 + 7A4 - As - 9Ag - A7 - Ag < 5 
-Ai + 7A2 - A3 - A4 + 7As - Ag - 9A7 - Ag < 5 
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-3Ai + 5A2 + 5A3 + I3A4 - IIA5 - SAg - IIA7 + SAg < 7 
5Ai + I3A2 - IIA3 + 5A4 - IIA5 - 3A6 - 3A7 + SAg < 7 
5Ai - 3A2 + 5A3 + I3A4 - IIA5 - llAe - 3A7 + SAg < 7 
5Ai + I3A2 - IIA3 - 3A4 + 5A5 - llAe - 3A7 + 5Ag < 7 

19Ai + IIA2 - 2IA3 - I3A4 - 5A5 - 5A6 + 3A7 + llAg < 9 
19Ai - I3A2 - 5A3 - 5A4 + 3A5 + llAe - 2IA7 + llAg < 9 
llAi + I9A2 - 2IA3 - 5A4 - I3A5 - SAe + 3A7 + llAg < 9 
-5Ai + 3A2 + IIA3 + I9A4 - 2IA5 - 13A6 - 5A7 + llAg < 9 

5Ai + A2 + A3 - 3A4 + A5 - 3A6 - 3A7 + Ag < 4 
Ai + A2 + 5A3 - 3A4 + As + Ae - 3A7 - 3Ag < 4 
Ai + A2 + A3 + A4 + 5A5 - 3A6 - 3A7 - 3Ag < 4 
Ai + 5A2 + A3 - 3A4 + As - 3A6 + A7 - 3Ag < 4 
5Ai - 3A2 + A3 + A4 + As + Ae - 3A7 - 3Ag < 4 
5Ai + A2 + A3 - 3A4 - 3A5 + Afi + A7 - 3Ag < 4 
5Ai + A2 - 3A3 + A4 + As - 3A6 + A7 - 3Ag < 4 

-Ai + 3A2 + 3A3 - A4 + 3A5 - Ae - A7 - 5Ag < 4 
3Ai + 3A2 - A3 - A4 + 3A5 - 5Ae - A7 - Ag < 4 
3Ai + 3A2 + 3A3 - 5A4 - As - Ae - A7 - Ag < 4 
3Ai - A2 + 3A3 - A4 + 3As - Ae - 5A7 - Ag < 4 
3Ai + 3A2 - A3 - A4 - As - Ae + 3A7 - 5Ag < 4 
3Ai - A2 - A3 + 3A4 + 3As - Ae - A7 - 5Ag < 4 
3Ai - A2 + 3A3 - A4 - As + 3Ae - A7 - 5Ag < 4 

4-3.6 Remark. The marginal inequahties are independent and written in the form ()4.3|1 
of theorem 14. 2. II Using the normahzation equation Tr p = n they can be transformed in 
many different ways. For example, the above constraints for system A^Tij are equivalent 
to inequahties ()4.2|1 . The inequalities for A'^Tig can be recast into a nice form found 
experimentally by Borland and Dennis 

+ 1x21 + |x3| + \x4\ + Ixsl + |xe| + \xj\ < 4, 

where 

Xi = Ai + A2 + A3 + A4 — As — Ae — A7 — Ag 
X2 = Ai + A2 + As + Ae — A3 — A4 — A7 — Ag 
3^3 = Ai + A3 + As + A7 — A2 — A4 — Ae — Ag 
X4 = Ai + A4 + Ae + A7 — A2 — A3 — As — Ag 
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X5 


= A2 


+ A3 


+ Ae 


+ A7 


-Ai 


-A4 


-As 


-A 




= A2 


+ A4 


+ As 


+ A7 


-Ai 


-A3 


-As 


-A 


X-j 


= A3 


+ A4 


+ As 


+ Ae 


-Ai 


-A2 


-A7 


-A 



Borland and Dennis numerical data were inconclusive for the system A^Tig described 
by 31 inequalities. One may wonder whether they can be written in a nice symmetric 
form. 
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